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Strangeness production in heavy ion collisions 

^ ! Krzysztof Redlich a * t 

a Gesellschaft fiir Schwerionenforschung, D-64291 Darmstadt, Germany 

Strangeness production in heavy ion collisions is discussed in a broad energy range 
from SIS to RHIC. In the whole energy range particle yields are showing high level 
of chemical equilibration which can be described by the unified freezeout conditions of 
fixed energy/particle ~ lGeV. The statistical model within the canonical formulation of 

■^j- ! strangeness conservation provides a framework to describe the observed enhancement of 
(multi) strange particles from p+A to A+A collisions measured at the SPS energy and 

iy-j ' predicts that this enhancement should be larger for decreasing collision energy. However, 
only at the SPS and RHIC chemical freezeout temperature is consistent within error with 
the critical value required for deconfinement and simultaneously strangeness is uncorre- 
! lated and distributed in the whole volume of the fireball. 

a: 

Ph! 1. Introduction 

<D . 

Ultrarelativistic heavy ion collisions provide a unique opportunity to study the prop- 
.£h ! erties of highly excited hadronic matter under extreme conditions of high density and 
^ ■ high temperature |l]-[5| . From the analysis of rapidity distribution of protons and of their 
transverse energy measured in 158 A GeV/c Pb+Pb collisions an estimate of the initial 
conditions @-|5[ leads to an energy density of 2-3 GeV/fm 3 and a baryon density of the 
order of 0.7/fm 3 . Lattice QCD at vanishing baryon density suggests that the phase tran- 
sition from confined to the quark-gluon plasma (QGP) phase appears at the temperature 
T c = 173 ±8 MeV which corresponds to the critical energy density e c ~ 0.6 ±0.3 GeV/fm 3 
||. One could thus conclude that the required initial conditions for quark deconfinement 
are already reached in heavy ion collisions at the SPS energy. It will be discussed to what 
extent the composition of the final-state hadrons, in particular their strangeness content 
can be considered as a probe of quark deconfinement in the initial state. 



2. Strangeness content of the quark-gluon plasma - signal for deconfinement 

An enhanced production of strange particles was long suggested as a possible signal of 
the QGP formation in heavy ion collisions [[j]|§. In the QGP the production and equilibra- 
tion of strangeness is very efficient due to a large gluon density and a low energy threshold 
for dominant QCD processes of ss production. In hadronic systems the higher thresh- 
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Figure 1. Ratio of total number of kaons 
per pions versus the number of participant 
in p+p, p+A and A+A collisions [pi] . 
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Figure 2. The multiplicity of K + , K and 
p per pion multiplicity function of the 
number of collisions in p+A reactions H. 



old for strangeness production was argued to make the strangeness yield considerably 
smaller and the equilibration time much longer. [] 

On the basis of the above strangeness QGP characteristics the following experimental 
implications are commonly quoted as the signal of deconfmement: 

i) global strangeness enhancement: strangeness content of secondaries, measured by the 
total number of produced < ss > quarks per participant A part or per produced light 
quarks < uu + dd> should increase from pp, pA to AA collisions. 

ii) enhancement of multistrange baryons: the specific enhancement (increasing with 
strangeness content) of multi-strange baryons and anti-baryons in central AA collisions, 
with respect to proton induced reactions follows deconfmement. 

iii) chemical equilibration of secondaries: the appearance of the QGP being close to 
chemical equilibrium and subsequent phase transition should in general drive hadronic 
constituents produced from hadronizing QGP towards chemical equilibrium. 

Heavy ion experiments at CERN SPS reported actually a global strangeness enhance- 
ment from p-p, p-A to A-A collisions fl|,|9|-|TT|j . In Fig. 1 we show full phase-space data of 
NA35 and NA49 collaboration on<K + K>/<n> ratio in S+S, S+Ag and Pb+Pb 
relative to p+p and p+A collisions. There is indeed a jump of factor of about two in 
strangeness content when going from p+p to heavy ion collisions. However, strangeness 
enhancement is already present in p+p and in p+A collisions |§. In Fig. 2 the 
< K + > / < 7r + > ratio measured by the NA49 collaboration in p+Pb collisions in 
the forward hemisphere shows an increase with the number of collisions due to the sec- 
ondary production. This effect is qualitatively similar to the one observed in Pb+Pb 
collisions. Strangeness enhancement from p+p to the most central p+A collisions can 
already account for more than 50% increase of strangeness seen in Fig. 1. 

Large strangeness content of the QGP plasma should be reflected in a very specific 
hierarchy of multistrange baryons J7|: enhancement of Q > H > A. Fig. 3 shows the 
yield/participant in Pb+Pb relative to p+Pb collisions measured by the WA97 and NA57 

3 In [fj9| it was argued, however, that multi-mesonic reactions could accelerate the equilibration time of 
strange antibaryons especially when the hadronic system is hot and very dense. 
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Figure 3. Particle yields per participant in Pb+Pb relative to p+Be and p+Pb collisions 



centrality dependence. The data are from WA97 [11[] and NA57 [[L2|] collaboration. 



collaboration [ TTj , 12 1 . Indeed the enhancement pattern of the (anti)hyperon yields is 
seen to increase with strangeness content of the particle and there is a saturation of this 
enhancement for N woni i > 100. Recent results of the NA57 collaboration are showing in 
addition an abrupt change of anti-cascade enhancement for a lower centrality. Similar 
behavior was previously seen on the level of K + yield measured by the NA52 experiment 
in Pb-Pb collisions ]13|]. These results are very interesting as they might be interpreted 
as an indication of the onset of a new dynamics. However, a more detailed experimental 
study and theoretical understanding are still required here. It is e.g. not clear what is an 
origin of different centrality dependence of 5 and S as well as inconsistences between the 
NA52 and the WA97 data. 
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Figure 4. The experimental data for particle 
multiplicity ratios at the SPS (full points). 
The predictions of thermal model (horizontal 
lines) 0. 



Figure 5. =S P article multiplicitysratios from 
STAR collaboration at RHIC (full points) 
f23| . The thermal model results are shown 



by the horizontal lines. 



A number of different mechanisms was considered to describe the magnitude of the 
enhancement and centrality dependence of (multi) strange baryons measured by the WA97 
[T4| |Tgf| . Microscopic transport models make it clear that the data shown in Fig. 3 can 
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not be explained by pure final state hadronic interactions. Only the combination of the 
formers with an additional pre-hadronic mechanisms like baryon junction processes, color 
ropes or color flux tubes overlap can partly explain the enhancement pattern and the 
magnitude for the most central collisions. However, the detailed centrality dependence 
is still not well reproduced within the microscopic models. An alternative description of 
multistrange particle production was formulated in terms of macroscopic thermal models 
and will be discussed in the next sections. 

3. Equilibration and particle yields 

One of the key questions in the description of particle production in heavy ion collisions 
is to what extent measured yields are showing equilibration. The level of equilibration 
of secondaries in heavy ion collisions can be tested on a different level: by analyzing (i) 
particle abundances or (ii) particle momentum distributions. In the first case one estab- 
lishes the chemical composition of the system (chemical freezeout), while in the second 
an additional information on dynamical evolution and collective flow can be extracted 
(thermal freezeout). 




Figure 6. Particle multiplicities per partici- 
pant normalized to its value in p+p system 
as a function of A vart . calculated in statistical 
model in (C) ensemble f21 



Figure 7. K + /tc + ratio at midrapidity in 
A+A relative to p+p collisions. For the 
compilation of data see |3"3" , |54]| . The full 
line is the parameterization from [33 



The compilation of the experimental data for particle multiplicity ratios in central 
Pb+Pb collisions at the SPS is shown in Fig. 4. The relative hadron abundances are 
compared with the thermal model [20]. The model was formulated in the grand canoni- 
cal (GC) ensemble with the partition function which contains the contributions of most 
hadrons and resonances and preserves the baryon number, strangeness and charge con- 
servation. The particle ratios depend only on two parameters; temperature T and baryon 
chemical potential \xb- It is clear from Fig. 4 that with T ~ 170MeV, corresponding to the 
energy density ej ~ 0.6GeV//m 3 , and with the baryon chemical potential /xg ~ 270MeV 
the statistical model reproduces the experimental data. 
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In Fig. 5 we compare the recent data of STAR collaboration in Au+Au collisions at 
RHIC with the statistical model results. Also here the model is consistent with data 
p2 , 24 1 -n The crucial difference to the SPS is a much lower value of the baryon chemical 
potential fiB of 45 MeV. This is to be expected as higher collision energy should lead 
to less stopping. The freezeout temperatures in central A+A collisions at the SPS and 
RHIC coincide within errors with the critical temperature from lattice QCD. This could 
indicate that all particles are originating from deconfmed medium and that the chemical 
composition of the system is established during hadronization 



Pb+Pb(Si-A) P.BraunMunzinger, et al 

Pb+Pb M.Gorenstein, et al 

Au-Au, Pb+Pb F.Becattini, A. Keranen, et al. 

SIS J. Cleymans, et al 

RHIC Nu Xu, et al 
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Figure 8. Average temperature and 
baryon chemical potential required to 
describe particle yields measured at 
GSI/SIS, AGS, and SPS [|5],g6|. The uni- 
fied chemical freezeout curve (dotted line) 
is calculated from the condition that the 
ratio of the average energy < E > per av- 
erage number of hadrons < N > is equal 
to 1 or 1.1 GeV. 



Chemical equilibration of secondaries is, however, not a unique property of the SPS 
data [pl],pE|-27|. It is also present at much lower incident energy or in peripheral high 
energy nucleus- nucleus collisions where the initial conditions required for deconfinement 
are not necessarily established. To discuss equilibration of strangeness in the above cases 
one needs, however, to formulate a statistical model in the canonical (C) ensemble with 



respect to the strangeness conservation [21,26,28,29 



4. Chemical equilibration - canonical description 

Strangeness conservation in statistical models can be described in the (GC) ensemble 
only if the number of produced strange particles per event is much larger than one. In the 
opposite limit of rare particle production pO| , strangeness conservation must be imple- 
mented locally on an event-by-event basis i.e., (C) ensemble of strangeness conservation 
must be used. The (C) ensemble is relevant in the statistical description of particle pro- 
duction in low energy heavy ion ||25|| , or high energy hadron-hadron or e + e~ reactions 
[p9"| as well as in peripheral heavy ion collisions [^TJ. The exact conservation of quantum 
numbers, that is the canonical approach, is known to reduce severely the phase-space 
available for particle productions. . 

To illustrate the canonical suppression let us consider a simple example of K + K~ 
production in the pionic thermal system of the volume V at temperature T. For large 
T, kaons are abundantly produced and the density of kaons reaches (GC) equilibrium 



4 A precise determination of Tf is not yet possible as fits with 170 < Tf < 190 give similar value of x 2 - 
Data like eg. A/p are necessary to further determine Tf [ pi) . 
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Figure 9. Statistical model results on the av- 
erage number of kaons Mr per participant 
Apart pfifl - Open (full) symbols represent 
measured multiplicities. 



Figure 10. K + / K data for different cen- 
trality at GSI/SIS H|, AGS B SPS 1 
and RHIC P5| . The dotted lines are sta- 
tistical model predictions 



value, = (l/27r 2 )m 2 K+ TK 2 (m K + /T) . In the limit of low temperature the K + ,K~ are 
very rarely produced and in order to satisfy strangeness neutrality condition, kaons must 
appear in pairs in the near vicinity. In this case, the density of produced kaons reaches 
an equilibrium result of (C) ensemble. In the asymptotic limit with < K ><< 1 we have 

n c K+ ~ [(l/27c 2 )m 2 K+ TK 2 (m K+ /T)] x [Vo(l/2n 2 )m 2 K .TK 2 (m K - /T)]. (1) 

The first term coincides with (GC) value, the second describes the phase space suppression 
since with each K + a K~ has to appear in the near vicinity in order to conserve strangeness 
exactly. The parameter Vq in Eq. 1 is introduced as a correlation volume where the K + 
and K~ should be created to fulfill locality of strangeness conservation. In heavy ion 
collisions V was found to scale with A part and in proton induced precesses V ~ V proton . 
There are thus, two origins of canonical suppression of strangeness: first, particles are 
produced in pairs what restricts the available momentum phase-space and secondly they 
appear in the near vicinity in space to fulfill locality of the conservation laws. 

4.1. Multistrange baryons - centrality dependence 

The canonical suppression of thermal particle phase-space increases with strangeness 
content of the particle. The exact conservation of strangeness requires that each particle 
carrying strangeness s has to appear e.g. with s other particles of the strangeness one to 
satisfy strangeness neutrality condition. 

In Fig. 6 we calculate the multiplicity/participant of Q, S, and A relative to its value in 
a small system with only two participants. Thermal parameters were assumed here to be 
Ap art independent. Fig. 6 shows that the statistical model in (C) ensemble reproduces the 
basic features of WA97 data: the enhancement pattern and enhancement saturation for 
large A part indicating here that (GC) limit is reached. The quantitative comparison of the 
model with experimental data would require an additional assumption on the variation 
of He with centrality to account for larger value of B/B ratios in p+A than in Pb+Pb 
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Figure 11. The K + /tc + , Att ratio versus en- 
ergy |34] . The short-dashed and dashed lines 
represent thermal model result for Pb+Pb 
and p+p respectively. The full lines are a 



parameterization of p+p data [33] 
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Figure 12. As in Fig. 11 but at midra- 
pidity. The results at RHIC were esti- 
mated from STAR results on K~ and 



K + 1 K ratios |23], assuming that the ra- 
tio Of 7T + /tT~ ~ 1. 



collisions | IT[| . An abrupt change of the enhancement seen in the NA57 data for 5 is, 



however, very unlikely to be reproducible in terms of this approach. 

One of the consequences of the model is that the enhancement pattern seen in Fig. 
6 should not be only a feature of the SPS data. In terms of (C) model strangeness 
enhancement and enhancement pattern should be also present there in heavy ion collisions 
at lower energies and should be even more pronounced. This is in contrast e.g. to UrQMD 



predictions which are showing increasing enhancement with beam energy |17 



In Fig. 7 we show a compilation of the data on K + /n + ratio in A+A relative to 
p+p collisions. This double ratio could be referred to as a strangeness enhancement 
factor. The enhancement is seen to be the largest at the smallest beam energy and 
is smoothly decreasing towards higher energy. If strangeness enhancement is indeed of 
thermal origin then similar behavior is also expected for multistrange baryons. This could 
put in question |JT the observed strangeness enhancement from p+p to A+A collisions 
as an appropriate characteristics of the deconfinement. To search for the QGP formation 
through the strangeness composition of secondaries one should rather look for a non- 
monotonic behavior of strange particle yields versus centrality or collision energy. 

A detailed analysis of the experimental data in heavy ion collisions from SIS to SPS 
pO| , |25| -F2^1 has shown that most particle yields are reproduced by the statistical model.Q In 
Fig. 8 we present the compilation of chemical freezeout parameters required to reproduce 
measured particle yields at SIS, AGS, SPS and RHIC energies. The GSI/SIS results have 
the lowest freezeout temperature and the highest baryon chemical potential. As the beam 
energy is increased a clear shift towards higher T and lower fiB occurs. There is a common 
feature of all these points namely that the average energy per hadron is approximately 1 
GeV. Chemical freezeout in A+A collisions is thus reached when the energy per particle 
drops below 1 GeV aX all collision energies [25]. Fig's. 9-12 are showing different particle 
multiplicity ratios along unified freezeout curve in comparison with experimental data. 



Statistical models, however, failed e.g. to describe the yield of r] at SIS or A/p ratio at AGS. 
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Figure 13. Dependence of the transverse Figure 14. Thermal ||,|37| and chemi- 
mass inverse slope parameters on the particle cal freezeout conditions p5| , p6| compared 
rest mass |]32| , ^0| . with the evolution path from UrQMD [jig] . 



In all cases the energy dependence of data is well reproduced by the statistical model. 
The canonical strangeness suppression in central A+A collisions is very relevant at SIS 
energies. At the top AGS the (GC) treatment of strangeness conservation is already 
adequate. 

Specific features of K + and K~ excitation functions shown in Fig. 9 are appearing in the 
model as a consequence of different dependence of K + and K~ yields on baryon number 
density ||26|| . Also here, the increase of the canonical correlation volume from the size of 
the colliding nucleus at SIS up to the macroscopic volume of a thermal fireball at SPS is of 
crucial importance. Furthermore, the independence of the K + / K~ ratio on the number 
of participating nucleons seen in Fig. 10 from SIS to RHIC is a natural result of the 
statistical model as a particle and its antiparticle have the same centrality dependence as 
well as in (GC) and in (C) formulation. A consistent description of the energy or centrality 
dependence of kaon yields in A+A collisions within microscopic transport models is still 
missing f35fl . 

Of particular interest are experimental results represented in Fig. 11 where the 4/T data 
on K + /tt + ratios |J4| are showing a maximum at ^/s ~ 6 — 8GeV and drop at the 
SPS. The suppression of this ratio disappears, however, at midrapidity |54| as seen in 
Fig. 12. Statistical model favors the behavior seen in Fig. 12. The description of fully 
integrated data at the top SPS energy is only possible in terms of the statistical model 
when introducing an extra parameter which accounts for strangeness undersaturation [^7 . 
The suppression of K + /n + ratio seen in Fig. 11 was predicted as a possible signal for 
deconfinement |BH]. 



5. From chemical to thermal freezeout 

Transverse mass distribution of strange and non strange particles clearly indicates in 
addition to the thermal one also the collective transverse flow component [jj],[5l,|37|l . In 
Fig. 13 the resulting increase of the transverse mass slope parameter T app with the rest 



mass of the particle is typical for transversely expanding source. Deviations of f2 [^ and 
J/^ seen in Fig. 13 most likely reflect their earlier decoupling from the system due to a 
small rescattering cross section with the surrounding medium. 

The compilation of thermal freezeout parameters from SIS to SPS, extracted from 
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particle spectra and two particle momentum correlation are shown in Fig. 13. The 
system seems to expand transversely from the chemical to thermal freezeout following 
the path of constant entropy (S) per baryon (B) [^]. The abundances of secondaries 
during this evolution are frozen in. In the statistical approach this can only happen if 
the phase space of particles like e.g. pion , kaon and nucleon are over-saturated by 
an effective chemical potentials with typical values of 60-120MeV at T t h ~ HOMeV . The 
evolution path and the moment of dynamical decoupling of the fireball is seen in Fig. 14 
to coincide with UrQMD description fll4,15|. The chemical freezeout curve of constant 



< E > / < N >~ lGeV appears in UrQMD |TJj when the system reaches the kinetic 
equilibrium (full squares in Fig. 14). 



6. Conclusions 

One of the most puzzling results in heavy ion collisions is presumably an observation 
that particles seem to be produced according to the principle of maximal entropy. In a very 
broad energy range from ^Js ~ 2.5 up to 130 GeV hadronic yields and their ratios observed 
in heavy ion collisions resemble chemical equilibrium population along unified freezeout 
curve determined by the conditions of fixed energy/hadron ~ lGeV. Strangeness follows 
this systematics and since their production mechanism is quite different at low and at high 
energies it seems that particle yields are losing information about dynamics. However, 
there are some characteristic features of the system at chemical freezeout in central A+A 
collisions being only present at SPS and RHIC: i) chemical freezeout temperature is 
within errors consistent with the critical temperature and ii) strangeness is un-correlated 
and redistributed in the whole volume of a thermal fireball. Although the condition i) 
is also there in peripheral nucleus-nucleus and even in p+p collisions, the property ii) is 
obviously not valid any more. Here, the thermal phase space available for strange particles 
is strongly suppressed since, for only a few particles being produced per event, strangeness 
is strongly correlated in the volume Vo ~ V pro ton and thus has to be conserved exactly 
and locally (canonical description). In high energy A+A reactions, there are already 
sufficiently many strange particles being produced per event such that strangeness is 
conserved on the average and is distributed in the whole volume of the fireball (grand 
canonical description) . 

If thermal description of particle productions is correct, strangeness enhancement from 
p+p to A+A collisions appears naturally as a transition from the canonical to the grand 
canonical regime. Hence strangeness enhancement does not necessarily require deconfine- 
ment and should also be observable at lower collision energies where the conditions for de- 
confinement are not satisfied. However, a specific, non-monotonic behavior of strangeness 
versus centrality or collision energy could be of interest in this context. Recent results of 
NA57 and NA52 showing an abrupt change with centrality of S and K + , respectively, as 
well as the observed by NA49 drop of K + /tt + ratio seen in Fig. 11 could signal a new 
dynamics. 
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